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' ' The energy region spanning from '-^ 10^^ to < 10^^ eV is critical for understanding 

^H ■ both, the Galactic and the extragalactic cosmic ray fluxes. This is the region where 

c/3 , the propagation regime of nuclei inside the Galactic magnetic environment changes 

from diffusive to ballistic, as well as the region where, very likely, the most powerful 

Galactic accelerators reach their maximum output energies. In this work, a diffusion 

^ . Galactic model is used to analyze the end of the Galactic cosmic ray spectrum 

(^ [ and its mixing with the extragalactic cosmic ray flux. In particular, we study the 

conditions that must be met, from the spectral and composition points of view, by 

.^ . the Galactic and the extragalactic fluxes in order to reproduce simultaneously the 

I I I total spectrum and elongation rate measured over the transition region by HiRes 

f^ ' and Auger. Our analysis favors a mixed extragalactic spectrum in combination with 

0|0 , a Galactic spectrum enhanced by additional high energy components, i.e., extending 

beyond the maximum energies expected from regular supernova remnants. The two 

additional components have mixed composition, with the lowest energy one heavier 

^ ■ than the highest energy one. The potential impact on the astrophysical analysis of 

H , the assumed hadronic interaction model is also assessed in detail. 
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Introduction 



The cosmic ray (CR) energy spectrum extends for many orders of magnitude 
with a power law index ^2.7. Along this range of energies, the three spec- 
tral features are known: the first knee at -E ~ 3 PeV, the second knee at 
E ft! 0.5 EeV and the ankle, a dip extending from the second knee to beyond 
10 EeV. The nature of the second knee and of the ankle is still uncertain 
[l]; a possible interpretation of the two features is the transition between the 
Galactic and extragalactic components. At energies between 10^^ — 10^^ eV 
the Galactic supernova remnants (SNR) are expected to become inefficient as 
accelerators. This fact, combined with magnetic deconfinement, should mark 
the end of the Galactic component of cosmic rays, although the picture could 
be confused by the existence of additional Galactic accelerators at higher en- 
ergy. On the other hand, at energies above the second knee, extragalactic 
particles are able to travel from the nearest extragalactic sources in less than 
a Hubble time. Consequently, the spectrum may present above 10^^'^ eV a 
growing extragalactic component that becomes dominant above 10^^ eV. The 
region between the second knee and the ankle could be the transition region 
between the Galactic and extragalactic components. 



In this work we analyze the transition region comparing the diffusive Galactic 
spectrum from SNRs with two different models of extragalactic spectrum, one 
in which only protons [2] are injected at the sources and another in which 
a mixed composition containing heavy nuclei [2| is injected. To discriminate 
between the possible astrophysical scenarios we analyze the composition pa- 
rameter Xjnax inferred from the different extragalactic (EG) models using the 
hadronic interaction models used in literature. A comparison with the exper- 
imental data is done with special attention to Auger and HiRes data. 



1 Galactic and extragalactic cosmic rays: transition models 



The nature of the knee, consisting of a steepening of the spectrum from a 
power low index jg = 2.7 to 7^ = 3.1, is well explained by the rigidity- 
models (rigidity-acceleration model and rigidity-confinement model), in which 
the maximum energy achievable by nuclei is rigidity dependent. In these mod- 
els, the knees of the spectrum of nuclei of charge Z are related to the proton 
knee energy E^^ = ZE^^, where E^^ = 2.5 x 10^^ eV. Beyond the highest en- 
ergy knee E[^ = 6.5 x lO^^eV^ the total Galactic fiux, dominated by the Iron 
component, must be steeper. While this spectral feature seems to be confirmed 
by KASCADE data [3], the nature of the second knee and of the ankle are still 
uncertain [T]. Although the region between the second knee and the ankle is 



naturally considered as the transition between the Galactic and extragalactic 
component, the way the transition takes place is very model dependent. The 
main models describing the transition are the ankle model, the dip model and 
the mixed composition model. 

In the ankle model the transition between the Galactic cosmic rays (GCR) 
and extragalactic cosmic rays (EGCR) occurs at the ankle i?a ~ 1 x 10^^ eV 
[5]|6H7.8.9yl0illfl2] . The dip in the CR spectrum is explained as the intersec- 
tion of a flat extragalactic component with the steep Galactic component. 
An advantage of this model is the extragalactic fiat generation spectrum E'"^ 
which provides reasonable luminosities of the sources. On the other hand, the 
ankle model doesn't work well in the framework of the rigidity-model. If the 
transition from Galactic to extragalactic CRs occurs at the ankle, other ad- 
ditional Galactic mechanisms are required to accelerate particles at energies 
beyond E^^ = 6.5 x lO^^el/. 

A different explication of the dip and a lower transition energy comes from 
the dip model [T3ll2] . In this model the dip is caused by e~^e~ pair production 
by the extragalactic protons with the CMB photons. The transition energy 
is determined by the energy at which adiabatic energy losses are equal to 
pair production energy losses. In this case the transition between the Galac- 
tic component and the extragalactic component takes place at lower energies, 
around the second knee, in agreement with the rigidity-model. The calculated 
position and shape of the dip, confirmed by Akeno-AGASA [H], HiRes [15], 
Yakutsk p^pT] and Fly's Eye [H] experimental dataLiJ, is "universal" as it 
doesn't depend on the type of propagation, on the source density and separa- 
tion (for separation distance d < 100 Mpc) and on the evolutionary model of 
the sources. The only factor that affect the dip is a heavy component of the 
EG spectrum [20|2T] . The shape and position of the dip is in agreement with 
observations for an EG pure proton composition with a maximum contami- 
nation of Helium nuclei of the order of 10%. 

On the other hand, as the dip model requires a generation spectrum E~'^''^, a 
solution is needed to prevent the too high emissivity needed at lower energies 
[T31I221I21I231. 



An alternative and intermediate model is the mixed composition model [2^f25]|3] . 
In this model the EG cosmic ray composition is assumed to be mixed, in anal- 
ogy with the Galactic component. As in the ankle model, the intersection of 
the Galactic and extragalactic component gives origin to the the dip structure 
but with the advantage of a lower transition energy (around E ^ 3x 10^^ ^V), 
which softens the requirement of additional acceleration mechanism. The cal- 
culated spectrum, fitted with the observed data, presents a spectral index 



^ Auger data don't contradict the high energy part of dip but an extension of Auger 
to lower energies it's needed to confirm the feature of the dip [T9] . 



7 ~ 2.1 — 2.3, which also provides a reasonable luminosity. 

The three model can be experimentally distinguished through measurements of 
the spectrum and of anisotropies, although the most discriminant feature is the 
chemical composition. While in the ankle model the Galactic heavy component 
dominates up to the ankle energy, in the dip model the transition is completed 
at energy around 1 x 10^^ eV where the composition is proton dominated. 
The composition at this energy (Iron/proton) is a discriminant between the 
two models. In the case of the mixed model, the transition occurs at E ^ 
3 X 10^^ eV and the chemical composition in the dip region is mixed. This 
model predicts a slower decrease of the Fe component and a slower increase 
of the proton fraction in the transition energy range. 



2 Galactic-extragalactic spectrum 



The current paradigm for Galactic cosmic ray acceleration is the Fermi accel- 
eration mechanism by shock waves of SuperNova Remnants (SNRs) [26] . 
In this section we calculate the Galactic diffusive spectrum from SNRs using 
the numerical diffusive propagation code GALPROP [271128] . 
The calculated Galactic spectrum from SNRs is then combined with two dif- 
ferent models of EG spectrum, one in which only protons [2] are injected at 
the sources and another in which a mixed composition containing heavy nuclei 
[5] is injected instead. 

We analyze the transition region between Galactic and extragalactic com- 
ponents in the two different EG scenarios and compare the combined total 
spectrum with the available experimental data. 



2.1 Diffusion Galactic model 



We use the numerical diffusive propagation code GALPROP [27f28] to repro- 
duce the Galactic spectrum from SNRs. 

The diffusive model is axisymmetric. The propagation region is, in cylindrical 
coordinates, bounded by i? = i?^ = 30 kpc and z = Zh = ^ kpc, beyond which 
free escape is assumed. The diffusion equation is: 

d4> - - 9 1 1 

-^ = q{f,p) + V X (D.^Vt/j) - —ijnlj) --4,--^ (1) 

at op Tf Tr 



where ip{f,p,t) is the density per unit of total particle momentum, q{r,p) 
is the source term, D^x is the spatial diffusion coefficient, p = dp/dt is the 



momentum loss rate and Tf and Tr are the time scale of fragmentation and 
the time scale of radioactive decay respectively. 



The GALPROP code solves the diffusion equation for all cosmic-ray species 
starting from the heaviest nucleus and then proceeding to lighter nuclei us- 
ing the computed secondary source functions. The numerical solution of the 
transport equation is based on a Crank-Nicholson [29P7] implicit second-order 
scheme. 



The diffusion coefficient is taken as (3Dq{p/ p^Y , where p is the particle rigid- 
ity, Do is the diffusion coefficient at a reference rigidity pD and 5 = 0.6. The 
diffusion coefficient can be inferred from the abundances of light nuclei, pro- 
duced mainly through spallation of heavy elements, as Li, Be and B, which give 
an estimation of the time of residence of CRs in the galaxy of ~ 1.5 x 10^ yr 
[30] . We use Dq = 5.75 x lO^^cm^s"^ at the reference rigidity p^ = A GV . 



The assumption of a diffusion coefficient with an energy dependence E~^'^ is 
not universally accepted [31j. In fact, the turbulence in the interstellar medium 
seems to follow closely a Kolmogorov spectrum, which should lead to an energy 
dependence of E^/'^ for the diffusion coefficient. In this scenario, the difference 
between primary and secondary cosmic ray energy spectra could be explained 
by the trapping of primary cosmic rays in high density regions with Kol- 
mogorov turbulence near their acceleration sites, where secondary CR would 
be mainly produced, and the diffusion of the latter inside this region and the 
interstellar medium. However, since at high energy both models produce simi- 
lar spectra and composition profiles, we think our analysis is to a large extent 
independent of these assumptions. 



The assumed distribution of cosmic rays sources is the one reproducing the 
cosmic-ray distribution determined by the analysis of EGRET gamma-ray 
data, which has the same parameterization of the R-dependence as that used 
for SNRs 



f{R,z) = f{R)e^^{-\z\/z,,aie) (2) 



where z scale = 0.2 kpc and: 



fiR) 



19.3 i? < 4 kpc, 

21.9 4kpc<R<8 kpc, 

15.8 8 A;pc < i? < 10 kpc, 

18.3 10 kpc< R< 12 kpc, 

13.3 12 kpc<R< 15 A;pc, 

7.4 i? > 15 A;pc. 



This choice is due to the fact that the SNR distribution [32j produces CR dis- 
tribution distinctly different from the measurements. A solution to the appar- 
ent discrepancy between the radial gradient in the diffuse Galactic gamma-ray 
emissivity and the distribution of SNRs was proposed by Strong et al. [33]. 

According to shock acceleration models [3HI35] . the injection spectrum is a 
power law function in rigidity with a break at rigidity po, beyond which it 
falls exponentially with a rigidity scale pc'- 
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P < Po, 
P> Po, 



where po = 1-8 PV, Pc = 1-26 PV and a = 2.05, which is the case of strong 
shock waves (M >> 1) [36ll37j . 



Stable nuclei with Z < 26 are injected, with energy independent isotopic 
abundances derived from low energy CR measurements |28ll38ll39f40ll41ll42f43j . 

The interstellar hydrogen distribution, molecular, atomic and ionized (H2, HI, 
HII), is derived from radio HI and CO surveys in 9 Galactocentric rings and in- 
formation on the ionized component [B] . The distribution of molecular hydro- 
gen is derived indirectly from CO radio-emission and the assumption that the 
conversion factor H2 /CO is the same for the whole Galaxy [35]. The atomic 
hydrogen (HI) distribution is taken from pB], with a z-dependence calculated 
using two approximation at different Galactocentric distance R [ITjHH]- The 
ionized component HII is calculated using a cylindrically symmetrical model 



The Interstellar Radiation Field (ISRF) is calculated using emissivities based 

on stellar populations and dust emission [50f5T] . 

We consider a Galactic turbulent magnetic field with a mean value of the 



component perpendicular to the CR propagation given by 

5 = Soexp(-^)expf-^), (3) 



where Bq = 6 /iG, Rscaie = 10 kpc, Zscaie = 2 kpc and Rq = 8.5 kpc the Sun 
Galactocentric distance. 



2.2 Diffusive Galactic spectrum 



The calculated diffusive Galactic spectrum is shown in Fig{T] superimposed 
to several experimental data results from the air shower experiments (HiRes 
[52J, HiResMIA [53], Fly's Eye ^ISj, KASCADE [4j, Yakutsk |I6], Akeno 
[5455J, AGASA [H], Haverah Park ^5m7\ . Auger [13:58], BLANCA [52|, Ti- 
bet [SOIIM] , Mt. Norikura [52]) as well as from direct measurements (SOKOL 
[63], JACEE [611I65] . Grigorov [66], Proton Satellite [67], Runjob [68]). 

It can be seen that there is considerable dispersion over the whole energy inter- 
val, which highlights the inherent difficulty in CR spectral measurements over 
so many decades in energy. Furthermore, it can be seen that, below the first 
knee the dispersion is mainly limited to flux normalization, while the spectral 
indexes seem to be fairly consistent for different experiments. The situation 
worsens at higher energies until large discrepancies are apparent above 10^^ 
eV. However, even at these high energies, a careful renormalization in energy 
[69] seems to be able to reconcile the several experimental data sets, giving 
a rather clear and consistent picture for the ankle. Unfortunately, normaliza- 
tion flux differences among spectra still remain between major experiments 
after the previous energy scale shifting, and they represent a concern to be 
experimentally addressed. 

The diffusive Galactic spectrum ^g is normalized to match KASCADE data 
at ~ 3 X 10^ GeV [4j. At this energy the differential flux value of the all 
particle energy spectrum and its uncertainty are the same for QGSJet 01 
and Sibyll 2.1 based analysis of the KASCADE experiment. While, with this 
renormalization, our spectrum agrees with data of various experiments at lower 
energies (JACEE, SOKOL, Tibet, KASCADE, Haverah Park and Akeno), 
beyond the knee the diffusive spectrum presents a strong deflcit of flux (Fig{T]). 

Since at -E > 10^ GeV the composition is dominated by intermediate {Z : 
6 — 12, i.e the CNO group) and heavier {Z : 19 — 26) nuclei, we renormalize 
these components {^q'^'^ , ^q) by a factor of 2, which produces a good agree- 
ment with the experimental data (see Figj2]). 
The renormalization of the CNO group and of the heavy nuclei group by a fac- 
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Fig. 1. Diffusive Galactic spectrum, ^g- The different Z-grouped nuclei components 
are shown. 
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Fig. 2. Diffusive Galactic spectrum, ^c, with renormalization of the CNO group 
(Z : 6— 12), ^q'^^, and of the heavy component {Z : 19 — 26), ^q. Agreement with 
experimental data is obtained for A:' ~ 1 and k ~ 0.8 — 1.2. 



tor of 2 is acceptable since it corresponds to a renormalization of the injected 
abundances into the acceleration mechanism. The Galactic CR CNO nuclei 
are over abundant with respect to Iron by several orders of magnitude. There- 
fore, the Iron contribution to CNO CR flux resulting from spallation is small. 
Consequently, to a good approximation, the observed Galactic abundances 
of Iron and the CNO group at the Solar circle can be varied independently 
by modifying their relative abundances at injection. This amount of relative 
renormalization is well inside the present uncertainties regarding the detailed 
workings of the injection mechanism. 

The renormalized diffusive Galactic spectrum reproduces well the data up to 
E pa 10*^ GeV, beyond which the spectrum falls steeply because of the end of 
the SNR acceleration efficiency. 



2.3 Extragalactic spectrum 



In order to study how the transition between the Galactic and extragalactic 
components takes place, we compare the Galactic spectrum originated in SNRs 
with two different possible scenarios for the extragalactic component. 

In the first model [2], a pure proton extragalactic spectrum, accelerated by 
a homogeneous distribution of cosmic sources, is considered. Local overden- 
sities/deficits of UHECR sources affect the shape of the GZK modulation, 
but do not affect the low energy region where the matching with the Galac- 
tic spectrum occurs. Within this model, we consider different cases of local 
overdensity/deficit of sources: 

(1) universal spectrum: n/uQ = 1; 

(2) overdensity of sources: n/rio = 2, R = 30 Mpc; 

(3) overdensity of sources: n/no = 3, R = 30 Mpc; 

(4) deficit of sources: n/no = 0, R = 10 Mpc; 

(5) deficit of sources: n/nQ = 0, R = 30 Mpc; 

where no is the mean extragalactic source density and n is the local overden- 
sity/deficit in regions of size R. 



In the second model |3j, the extragalactic spectrum is calculated for a mixed 
composition at injection typical of low energy cosmic rays. The spectral in- 
dex (3 is determined fitting the high energy CR data. Three different source 
evolution models in red shift are considered: 



a) strong evolution model: the injection rate is proportional to (1 + z)'^ for 
z < 1 and constant for 1 < z < 6, followed by a sharp cut-off, /5 = 2.1; 

b) SFR model: the EGCR injection power is proportional to the star formation 
rate which correspond to a redshift evolution (1 + zY for z < 1.3 and a 
constant injection rate for 1.3 < 2; < 6 (with a sharp cut-off a.t z = 6), 
13 = 2.2; 

c) uniform source distribution model: no evolution, (3 = 2.3. 



In both cases, the various parameters of the models are tuned to fit the avail- 
able CR data at UHE and are, in that highest energy regime, experimentally 
indistinguishable at present. The spectra of both the models, used in the next 
sections, are renormalized to HiRes data. 
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Fig. 3. Galactic and extragalactic spectrum matching for the proton models for an 
EG lower energy limit of 10^ GeV. The sum of the renormalized diffusive Galactic 
spectrum and of the extragalactic spectrum {^eg) for different cases of local over- 
density/deficit of sources ( §2.3p is shown for different renormalizations of the heavy 
component {^q)- The CNO group {^q^^) of the diffusive Galactic spectrum (<5g) 
has been renormalized by a factor 2.2. 

2.4 Combined spectrum: m,atching Galactic and extragalactic components 



In order to study how the transition between the Galactic and extragalac- 
tic components takes place, we subtract the combined theoretical (Galactic 
$G plus extragalactic $sg) spectrum from the available data. Two different 
approaches are used. 

First, we try to match the experimental data by varying the normalization of 
the heavy Galactic component $^, while keeping constant the previous renor- 
malization of the CNO group ^q^*^ ( §2.21) . The best reproduced spectrum for 
the two extragalactic models are shown in FigsjS], H] and [51 In the case of the 
proton model, a discontinuity appears when the two spectra are added, re- 
gardless of the lower limit adopted for the extragalactic component: 10® GeV 
(FigO]) or 5 X 10^ GeV (FigJl]). The latter corresponds to cosmic accelerators 
operating for the entire Hubble time. 

For both, proton and mixed- composition models, there is a flux deficit above 
10® GeV. The problem is much stronger for the mixed- composition model 
where, regardless of the luminosity evolution of EG CR sources, the total 
spectrum presents a large deficit of flux between 10® GeV and ?« 3 x 10^ GeV 

(FigED. 



In order to solve this flux deficit, the only way out seems to be the intro- 
duction of an additional Galactic component. We estimate this component by 
subtracting the sum of the diffusive Galactic and extragalactic fluxes from 



10 




Fig. 4. Idem to FiglS] but for the EG proton model with a lower energy limit of 
5 * 10^ GeV. 
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Fig. 5. Idem to FigHJbut for the Galactic and extragalactic spectrum for the mixed- 
composition model. The CNO group (<I>^^'^) of the diffusive Galactic spectrum ($g) 
has been renormalized by a factor 2. Three source evolution models are considered 
iUM- (a) strong, (b) SFR and (c) uniform. 



a smooth fit to the world data. This method confirms us the need of the 
renormalization of the CNO group by a factor ^ 2.2 and ~ 2 for the proton 
and mixed- composition models, respectively. In the case of the proton models 
(FigsinilT]), the observed deficit can be resolved with the introduction of an ad- 
ditional cosmic ray flux component (the first additional Galactic component, 
GAl) ^GAi- For the mixed- composition models this is not enough and one 
more additional cosmic ray flux component must be introduced (the second 
additional Galactic component, GA2) $ga2 (FigE]). The additional component 
^GAi common to both families of models is obtained with a shift in energy of a 
factor ?^ 3 of the diffusive Galactic heavy component, renormalized by a factor 
^ 0.8 in the case of the mixed- composition models and < 0.6 for the proton 
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Fig. 6. Galactic and extragalactic spectrum for the proton model for an EG lower 
limit of 10*^ GeV: the additional Galactic component (pcAi and the total high energy 
Galactic component ((/)g + 4>gai) are also shown. 
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Fig. 7. Idem to FiglH] but for the EG proton model with a lower energy limit of 
5 X 10^ GeV: the additional Galactic component (pcAi and the total high energy 
Galactic component (t/)^ + (J)gai) are also shown. 

models respectively. The second additional component $ga2 is obtained in an 
analogous way but with an energy-shift factor of ~ 10 and a renormalization 
by a factor ?» 0.2. The corresponding spectra are shown in FigsEl [7] and [HI 

2.5 Discussion on the spectrum 



We have analyzed the matching conditions of the Galactic and extragalactic 
components of cosmic rays along the second knee and the ankle. 
It seems clear that an acceptable matching of the Galactic and extragalactic 
fluxes can only be achieved if the Galaxy has additional accelerators, besides 
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Fig. 8. Galactic and extragalactic spectrum for the mixed- composition model: the 
two additional components (pcAi, 4>GA2 and the total high energy Galactic compo- 
nent {(I)q + 0GA1 + 4>GA2) are shown. 

the fiducial SNRs assumed here, operating in the interstellar medium. 

In the particular case of the proton model, only one additional component is 
required. This could well represent the contribution from compact and highly 
magnetized SNRs, like those occurring in the central, high density regions of 
the Galactic bulge, inside the dense cores of molecular clouds or those expand- 
ing into the circumstellar winds of their progenitors. Actually, this component 
does not need to originate in a particular type of source in itself, but could 
be the result of a non-homogeneous SNR population drawn from a spectrum 
of progenitor masses and evolving in different environments corresponding to 
the various gas phases that fill the interstellar medium [70]. From the point 
of view of CR luminosity this should not be a problem, since the Fe compo- 
nent and the additional high energy component amount to ~ 9% and < 1%, 
respectively, of the total Galactic diffusive spectrum. 

In the specific case of the proton model, it is apparent that a perfectly smooth 
match between the Galactic and the extragalactic spectra is not possible. 
Some discontinuity or wiggle seems unavoidable in the energy spectrum be- 
tween ~ 3 X 10^^ and ~ 10^"^ eV (see figures [3], HI [6] and [7]) . This is due to the 
fact that the extragalactic spectrum should have a rather abrupt low energy 
cut-off due to the finite distance to the nearest extragalactic sources and their 
limited age. Of course, whether such spectral feature is actually observable is 
strongly dependent in practice on the magnitude of the experimental errors in 
the determination of the primary energy from shower measurements and on 
the available statistics. 



The matching of the mixed- composition model, on the other hand, has wider 
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astrophysical implications for the Galaxy. The Galactic CR production has to 
be extended up to at least the middle of the dip and this requires, besides the 
previous additional component, another high energy Galactic component. The 
origin of these cosmic rays pushes even further the acceleration requirements 
imposed on the Galaxy. It is likely that a different population of Galactic 
accelerators must be invoked. Viable candidates could include, for example, 
rapidly spinning inductors, like pulsars or magnetars, or even very high energy 
episodic events like Galactic gamma ray bursts [7T]. If this were the case, it is 
very likely that photon emission at TeV energies should uncover the sources. 
Changes in propagation regime inside the Galaxy at these high energies should 
manifest as an anisotropic component embedded in the isotropic incoming 
extragalactic background. Again, considerable statistics might be necessary 
to make such effect observable. 

The energy requirements involved in the production of the second additional 
component are rather modest: ~ 10^^ of the SNR CR component, or ~ 10^^ 
erg/sec pumped into particles between ~ 10^^ and ~ 10^^ eV. Therefore, few, 
or even a single source, could be responsible for this component. This carries 
attached the potential problem of undesirable fine-tuning because of the re- 
quirement that, at this precise moment in time, the CR luminosity of these few 
(or this single) sources is such that it allows for a smooth spectral matching 
along the ankle. 



It is important to note that the region between some few times 10^^ eV and 
approximately 10^^ eV is a transition region from the point of view of propa- 
gation p!|72] . In fact, Larmor radius (in pc) of a CR nucleus of charge Z can 
be conveniently parametrized as: 

10 f EEev\ /.N 



where Eeev is its energy and B^g is some appropriate average of the Galactic 
magnetic field inside the propagation region. Since the transversal dimensions 
of the Galactic disk are on the order of a few times 10^ pc, we see that the 
diffusion approximation for protons should start to be broken somewhere be- 
tween the second knee and 1 EeV. The same should happen to heavier nuclei 
at progressively higher total energies: ~ 3 EeV for the CNO group and ~ 10 
EeV for Fe. The corresponding transition for each nuclei should be gradual, 
with the propagation eventually becoming ballistic at higher energies. This 
transition results in a complicated picture in which the end of the Galactic 
confinement spans almost 1.5 decade in energy, depending on the particular 
nuclei considered. Our interpretation of the additional Galactic components 
(but not their necessity) is affected to some extent by the implicit assumption 
that the particles propagate diffusively. This assumption is increasingly ques- 
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tionable for protons as we approach 1 EeV, but should be vahd for heavier 
nuclei in the energy range considered here. We are currently working on a 
detailed analysis of these effects which will be presented elsewhere. 



3 Composition 



The different extragalactic models are able to produce total spectra that are 
indistinguishable within the current experimental resolution. The composi- 
tion of UHECR is essential to understand the transition between Galactic 
and extragalactic cosmic rays and to discriminate the different astrophysical 
scenarios. 

Since UHECR experiments do not directly measure the composition, we have 
to infer it from parameters characterizing the shower development profile at 
ground. One of the most reliable parameters is X^ax, the atmospheric depth 
of the maximum longitudinal development of a shower [TS]. The variation of 
(Xmax) with energy gives information about the change in composition of the 
CR flux. 

Different hadronic interaction models (HIMs) can be used to interpret the 
Xmax dependence on energy and on primary composition, and this is the 
main uncertainty associated with this parameter. In this section we will use 
different parameterizations of (X^ax) deduced from the hadronic interaction 
models currently in use, in order to infer the composition energy profile for 
the different cases of Galactic-extragalactic combined spectrum described in 
section p^il 



3.1 Hadronic interaction models 



The composition predictions depend on the HIM used. For the sake of com- 
pleteness, we take into account in what follows all the HIMs currently used in 
the literature: 

• EPOS 1.6 

• QGSJet 01 

• QGSJetHv02 

• QGSJetHvOS 

• Sibyll 2.1 

For each hadronic interaction model, we calculate a parameterization of (Xmax) 
as a function of the hadronic interaction model and of the primary energy and 
composition: 
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Table 1 

Coefficients of the parameterization of {X^ax) as a function of primary energy and 
composition (eq. [5]) for different hadronic interaction models. 



Model 


a 


h 


c 


P 


a 


/? 


QGSJetII-v03 


-0.033 


-0.010 


58.313 


184.315 


9.60 X 10-^ 


-2.41 X 10-3 


QGSJetII-v02 


0.184 


-5.930 


115.111 


-15.895 


2.00 X 10-3 


6.50 X 10-2 


QGSJetOl 


0.074 


-3.695 


103.804 


-19.123 


6.85 X 10-^ 


2.14 X 10-2 


EPOS 1.6 


-0.011 


1.831 


29.212 


265.847 


-1.47 X 10"^ 


4.69 X 10-"^ 


Sibyll 2.1 


0.300 


-8.422 


134.960 


-72.051 


3.38 X 10"^ 


2.62 X 10-3 



{Xmax){E,A,l) 






1 + a,A) 



(5) 



+Piil + PiA), 



where E is the primary energy, A is the atomic mass of the primary CR, i is 
the hadronic interaction model, e ~ 81 MeV is the critical energy in air and 
Qi, bi, Ci, Pi, ai, /3i are coefficients determined by fitting equation [5] to shower 
simulations (see Tabled]). 

Parameterization [5]is a good approximation in the energy range [10^"^ eV, lO^^ eV] . 
Our parameterization for p and Fe, for several HIMs, is shown in Fig J9] together 
with the experimental results of HiRes [H], HiResMIA [53|, Fly's Eye [73] and 
Auger 



This parameterization allows us to compute the average X^ax for showers 
generated by any nucleus of interest and for all observed primary energies 
above the second knee. 

With this parameterization, in the following sections we reproduce (X^ax) 
energy profiles for different composition scenarios in the context of the various 
cases of the combined Galactic-extragalactic spectrum. 



3.2 Galactic-extragalactic combined spectra: X^ax energy profiles 



Heretofore, we have analyzed from the point of view of the energy spectrum, 
different scenarios for the combination of the Galactic flux with alternative 
extragalactic models. For each of these scenarios, and under different assump- 
tions for the composition of individual flux components, we estimate X^ax 
along the transition region. 
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Fig. 9. (Xmax) parameterization as a function of primary energy and composition. 
For each HIM, the {X^ax) dependence on energy for proton and Iron primary is 
shown. Experimental data from HiRes, HiResMIA, Fly's Eye and Auger are also 
shown. 



As seen in §2.41 an acceptable matching of the Galactic and extragalactic 
fluxes can only be achieved if the Galaxy has additional accelerators besides 
the regular SNRs assumed in §2.1[ This can be achieved by including either 
one (pure proton EG spectrum) or two {mixed- composition EG spectrum) 
additional Galactic components. 



In order to infer X^ax for the combined flux, we have to do some assumptions 
on the composition of the additional Galactic components. 

The first additional Galactic component (GAl), which is needed in both of 
the EG scenarios previously considered, is probably contributed by compact, 
highly magnetized SNRs. Consequently this component is likely to be domi- 
nated by heavy elements, say. Iron, at the highest energies. 

The second additional Galactic component (GA2), on the other hand, might 
be associated with a minor population of still more powerful SNR accelerators 
or with a completely different population of particle accelerators, e.g., rapidly 
spinning inductors, like pulsars or magnetars. In the first case, GA2 would 
be Fe dominated, while in the latter case, GA2 could be proton dominated. 
Consequently, we consider two different possibilities for GA2: pure proton and 
pure Fe. 
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Fig. 10. Mixed- composition EG scenario of Allard and co-workers. Contribution to 
the extragalactic flux of different nuclei in the "uniform source distribution model" . 
(Adapted from: [25j). 

We calculate X^ax as a function of energy in the energy range [10^''' eV, 10^° eV] 
for all the possible combination of Galactic and extragalactic spectra previ- 
ously discussed in §2.41 from the point of view of the shape of the energy 
spectruniil: 

• pure proton EG model: we consider the combination of the calculated Galac- 
tic spectrum from SNRs with the EG "universal" (see §2.31) spectrum in the 
two cases of lower EG energy limit (5 x 10^^ eVand ~ 10^"^ eV^); the compo- 
sition of the Galactic additional component is assumed to be heavy (Iron); 

• mixed EG composition model: we consider the combination of the calculated 
Galactic spectrum from SNRs with the mixed EG spectrum in its "uniform 
source distribution model" variant (see §2.31) : we also assume that GAl is 
composed by pure Fe, while we analyze two opposite scenarios for GA2: 
pure proton and pure Fe. The EG composition is taken from Fig{TO] [25] . 



^ For the diffusive Galactic spectrum from SNRs we kept the renormalization of 
the CNO group by a factor ~ 2 determined in the previous analysis ( ^2.4p 
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Fig. 11. Xmax, EPOS 1.6: the Xmax profiles calculated for all the possible combined 
spectra are shown superimposed onto the experimental data. The curves represent: 
(i) the combined spectra for an EG proton model with lower energy limit ~ 10^^ eV 
(cyan blue) and 5 x 10^^ eV (green) respectively; (ii) the combined spectra for an 
EG mixed- composition model with pure proton (purple) and pure Fe (brown-orange) 
GA2. (Xmax) values are calculated using the EPOS 1.6 HIM. 

For all these scenarios with different component spectra and composition as- 
sumptions, we calculate {Xmax) values in the energy range [10^^ eV, 10^° eV] 
using the parameterization given in §3.11 for the various HIMs in current use. 



The results are shown for each HIM in FigsJTTl [T2l [T3l [H] and [I5l together 
with the experimental results of HiRes [71], HiResMIA [53], Fly's Eye [75] and 
Auger 



First of all we notice that, as expected, for all the HIMs the theoretical X^ax 
values in the case of EG proton model are indistinguishable for the two dif- 
ferent EG lower energy limits (green and cyan blue curves are superimposed 
on the figures). On the other hand, the model is very much compatible with 
HiRes-MIA data below 10^^ eV while, at higher energies, (Xmax) is too large 
when compared to any of the available data sets. 

As expected, much more variability in the calculated elongation rates is ob- 
served for the lower energy regime in the case of EG mixed- composition models, 
depending on the assumptions made about the composition of GA2 (purple 
and brown-orange curves). However, at energies larger than ~ 3 x 10^*^ eV, 
all the solutions converge to the same profile since the effects of GA2 become 
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Fig. 12. Xmax, QGSJetOl: idem to FigHUbut using the QGSJetOl HIM. 
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Fig. 13. X^ax, QGSJetll v02: idem to FigHUbut using the QGSJetll v02 HIM. 

progressively negligible on the combined composition profile. All these results 
are qualitatively independent of the assumed HIM. Nevertheless, even if all 
the HIMs give basically the same trend at all energies, actually what data set 
is quantitatively compatible with the mixed composition theoretical models 
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Fig. 14. X^ax, QGSJetll v03: idem to FigHUbut using the QGSJetll v03 HIM. 
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Fig. 15. Xmax, Sibyll 2.1: idem to FigHUbut using the Sibyll 2.1 HIM. 

at energies beyond 3 x 10^^ eV does depend on the assumed HIM. Further- 
more, at the highest energies, > 10^^ eV, there could be an indication that 
the models give a systematically lighter composition than what the data sug- 
gests. Unfortunately, at present, statistics at these energies are too low for all 
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experiments as to render any solid conclusion. 

The most critical energy region for the understanding of the Galactic-extragalactic 
transition and disentangling its flux components, is ~ 10^^ - ~ 3 x 10^^ eV. 
Inside this region, different experimental results, and they are different indeed, 
imply very different astrophysical scenarios, in particular with regard to the 
highest accelerators present in our own Galaxy. 

In the following sections we will center on two of the most important experi- 
mental results today, those of Auger and HiRes, and test what, if any, further 
modifications should be applied to the composition profile in order to make 
the data and the theoretical models as compatible as possible. 



3.3 Composition evolution along the transition region 



If our understanding of (X^ax) in terms of composition is reasonably correct. 
Auger and HiRes elongation rate data suggest a composition profile along the 
ankle compatible with a mixed extragalactic composition, despite differences 
in the energy dependence of Xmax for both experiments. 

The agreement between Auger X^ax data [76] and our predictions is remark- 
able at energies around 3 x 10^^ eV for all HIMs. This agreement extends to 
higher energies, even beyond 2 x 10^^ eV, for QGSJetl and QGSJetll v02 and 
v03, while EPOS and SIBYLL display different trends. The picture is more 
complicated at lower energy where, under the present assumptions, there is 
no agreement regardless of the assumed HIM. 

The situation almost reverses in the case of the HiRes-MIA combined data 
set, for which a good agreement can be obtained at low energies with the EG 
proton model, but there is no clear fit at higher energies. 

It must be noted, however, that the composition profiles shown in FigsilT], [121 
[T3l UM and [151 result simply from the combination of our previous solution 
to the total energy spectrum for the mixed model and the unmodified mixed 
model compositions as determined by [2S], plus fixed compositions for GAl 
and GA2. 

Therefore, in this section we analyze whether the composition profile of the two 
additional Galactic components and of the extragalactic flux can be modified 
in a suitable way in order to satisfy, simultaneously, the existing spectral and 
elongation rate data along the transition region. 

We use to this end a simplified model in which the shape of the extragalactic 
energy spectrum is kept, to first order, identical to AUard's mixed composition 
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model [3], but the composition is limited to just two nuclei, proton and Iron, 
whose relative abundances can change appropriately in order to reproduce the 
behavior of (Xmax) as a function of energy. The admitted lack of consistency 
in this approach is, we believe, compensated by the insight gained into the 
phenomenological constraints imposed by the data on the astrophysical models 
at plajcJ. 

On the other hand, the spectral shapes of both additional components and 
their normalizations are preserved, but their compositions as a function of 
energy are now functions to be determined during the fitting process. A binary 
mixture of p and Fe is also assumed. 

We further assume, a priori, that the diffusive Galactic spectrum and its com- 
position are the ones determined previously in §2.41 with the calculated renor- 
malization of the CNO group. The diffusive Galactic flux from SNRs, spectrum 
as well as composition, is kept constant afterward during the fitting procedure. 

We apply this procedure to both. Auger and HiRes data in order to gain an 
insight on how present experimental uncertainties can affect our astrophysical 
understanding of the transition region. 



3.3.1 Auger data 

FigHSl shows the total spectrum and Galactic components that are used in 
this section. For each source evolution model (see §2.31) . we renormalize the 
EG spectrum to the surface detector (SD) Auger flux at 10^^°^ eV [58j, where 
one can expect a negligible contribution from the Galactic CR flux. After- 
ward, we renormalize the GAl and GA2 in order to match the total spectrum 
with the Auger spectrum data; the hybrid Auger spectruno at lower energy, 
10^^ — lO^'^'^ eV [To], is used to this end. The lowest energy branch of the 
ankle is approximated by a smooth interpolation between the highest energy 
KASCADE spectrum and the lowest energy hybrid Auger data. 

It must be noted that, regardless of the renormalization, the extragalactic 
theoretical spectral models are too soft to allow a good fit to the Auger data 
at the highest energies. This inconsistency have also effects at lower energies, 
where the position of the ankle is artificially pushed down from the position 
determined by Auger at ~ lO^*^'^ eV. Therefore, as a compromise, we consider 



^ An alternative analysis to the spectral fit is the use of slanted shower data which, 
for Akeno data was performed by |70j . 

^ The Hybrid technique consists in reconstructing the CR showers using the com- 
plementary information given by the two independent detectors, the fluorescence 
detector (FD) and the surface detector (SD). These data allow the reconstruction 
of the cosmic ray energy spectrum up to energies below the mid-ankle. 
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Fig. 16. Combined total spectrum ^Q^ + ^'^^f renormalized to Auger spectrum data 
for different EG evolution models ( §2.3p . The Galactic spectrum $g, the Galactic 
additional components ^gai, ^GA2 and the total Galactic spectrum <1>^* are also 
shown. The normalization of the second Galactic additional component is the one 
determined for the EG "uniform source distribution" model. See the text for details. 

here the model corresponding to a "uniform source distribution" , which is the 
hardest one and provides the best visual fit to the highest energy data. 

From the estimated combined Galactic-extragalactic total spectrum and using 
the Xmax parameterization given in §3.11 eq. 0, we compute the relative abun- 
dances of proton and Fe nuclei of the EG spectrum and of the two additional 
Galactic components (GAl and GA2) that match the Auger (X^ax) data. 
The EG, GAl and GA2 compositions are determined in the energy range of 



available Auger X^ 
interaction models. 



ciata, 



3 X 10^'' to ~ 5 X 10^^, for all the hadronic 



The corresponding Iron and proton fluxes for the EG spectrum and for GAl 
and GA2 are shown in FigJTTl The different curves for each nucleus correspond 
to inferences obtained with different HIMs. 



24 



«G +*GAI +«GA2 +*EG ' 



Auger Hybrid i — e — 
AjgerSD >-a-- 



QGSJelOl - 

QGSJetll v02 - 

QGSJetllu03 ■■ 

Brbyll2.1 ■■ 

EPOS 1.6 - 




■.-.•:-S' ' ^.i*i"iv:3*"' 



p 



EG mixed composition model: 
proton component 



Auger Hybrid >- 
Auger SD >■ 



QGSJetOI - 
QGSJetll v02 - 

QGSJetll v03 ■ 
Sil3yll2.1 ■ 

EPOS 1.6 ■■ 




T T . "...?.....■ 



#,.-■•■ 



Fe 



EG mixed composition model: 
Iron component 



Energy (eV) Energy (eV) 

(a) EG composition: proton (left) and Iron (right) components. 



™ 1 e+2' 






Auger Hybrid .— e — i 
Auger SD !■■&■■; 




i II II 

if-"- — 



GAl, GA2 mixed composition: 
proton component 



■^G + *'GA1 + 'l'GA2 + "I'EG ' 




i II II 

" a J.-';,' 

jl--- '■■"■"■" 



GAl, GA2 mixed composition: 
Iron component 



(b) GAl and GA2 composition: proton (left) and Iron (right) components. 

Fig. 17. Auger data, (a) EG proton and Iron composition; (b) GAl and GA2 proton 
and Iron compositions. The fluxes are calculated, for each HIM, in order to repro- 
duce the Auger (Xmax) data. The Galactic additional components {^gai, ^gA2), 
the EG spectrum {^eg) and the total spectrum {^tot = ^G + ^GAi + ^GA2 + ^eg) 
are also shown, as well as the Auger spectrum data. 

The corresponding average atomic weights, < A >, for each HIM are shown 
in Fig[ 



3.3.2 HiRes data 

In order to assess the astrophysical imphcations of the present experimental 
uncertainties in the total spectrum and elongation rate, we also apply here the 
same procedure as in the previous section to the HiRes data [71|52f53j . Fig. 
[T8] shows the total spectrum, for different assumptions regarding the redshift 
evolution of the sources, and the Galactic components under consideration. 
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Fig. 18. Combined total spectrum $^* + $^g^ renormalized to HiRes spectrum data 
for different EG evolution models (i ^2.3p . The Galactic spectrum <I>g, the Galactic 
additional components ^gai-, ^ga2 and the total Galactic spectrum $^* are also 
shown. The normalizations of the two Galactic additional component are the ones 
determined for the EG "SFR" model. 

Since the HiRes spectrum extends to lower energies than those of Auger, in 
this case we renormalize both additional Galactic components to match the 
HiResll data in the energy range lO^'^'-^ — 10^^-^ eV (Fig{TS]). The EG spectra, 
on the other hand, are already normalized to HiRes data ( §2.3p . 

As in the Auger case, the position of the ankle measured by HiRes ( 6 x 
10^*^ eV) is not accurately reproduced with any of the combined spectra if a 
reasonable agreement at higher energies is required. Consequently, we select 
for subsequent analysis the EG "SFR" model, which is the one that provides 
the best visual fit to the experimental data. 



We assume, as in §3.3.11 that the composition of the EG spectrum and of the 
two Galactic additional components is a binary mixture of proton and Iron 
nuclei and compute their relative abundances in order to match the {X^ax) 
energy dependence of the HiRes data. 
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(b) GAl and GA2 composition: proton (left) and Iron (right) components. 

Fig. 19. HiRes data, (a) EG proton and Iron composition; (b) GAl and GA2 
proton and Iron compositions. The fluxes are calculated for each HIM in or- 
der to reproduce the HiRes and HiResMIA {Xmax) data. The Galactic addi- 
tional components {^gai^ ^gA2)-, the EG spectrum ($eg) and the total spectrum 
{^tot = ^G + ^GAi + ^GA2 + ^eg) are also shown, as well as the HiRes spectrum 
data. 



The compositions as a function of energy of EG, GAl and GA2 are determined 
for all the HIM in the energy range of HiRes and HiResMIA X^ax data ~ 10^^ 
to ~ 2 X 10^9. 

The results for the calculated Iron and proton fluxes for GAl, GA2 and EG 
are shown in FiglT^ 

The average atomic weights, (A), are shown in FigJSOl for each HIM. 
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Fig. 20. Average atomic weight of the combined cosmic ray flux: the values were 
obtained from the calculated EG, GAl and GA2 compositions in order to reproduce 
the Auger (Fig [T7|) and HiRes (Fig llQ]) {Xmax) data, under different HIM assump- 
tions. 

3.4 Discussion on composition 



In section 132] we calculated the Xmax dependence on energy for all the possible 
scenarios obtained by combining the calculated SNR Galactic spectrum with 
different extragalactic models ( §2.4p . Several HIMs were used. 

Simple assumptions were made on the composition of the two additional 
Galactic components. The first one, which is required by both EG models, 
was assumed to be Iron dominated, as it is probably contributed by compact 
and highly magnetized SNRs. The origin and properties of GA2, required in 
the case of a mixed composition EG spectrum, are more uncertain. Conse- 
quently, as an exploratory test, we considered in section 13.21 two cases for 
GA2: a pure proton and a pure Iron composition. 

For all the possible combined spectra, we used different parameterizations of 
{Xmax) obtained for the HIMs used in literature to calculate the behavior of 
Xmax in the energy range [10^''' eV, 10^° eV]. We compared the calculated Xmax 
energy profiles with the available experimental data from Fly's Eye, HiRes, 
HiRes-MIA and Auger. 

The EG proton model is very much compatible with HiRes-MIA data at en- 
ergies below 10^^ eV for any HIM. However, if the current interpretation of 
Xmax in terms of shower composition is correct, the later model is completely 
incompatible with higher energy data which points to a mixed composition. 

In the case of the EG mixed-composition model, also two energy regimes are 
clearly distinguishable. 

At energies above ~ 3 x 10^*^ the slope of the elongation rate seems to be consis- 
tent with that of the data up to, possibly, ~ 2 x 10^^ eV. The effect of different 
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HIMs is basically to change the normalization of the curves, changing at the 
same time the experimental data set with which the models are more compati- 
ble. For energies smaller than ~ 3 x 10^^ eV, the calculated {X^ax) profiles are 
very dependent on the assumption made regarding the composition of GA2. 
A heavy GA2 can fit reasonably well the Fly's Eye data, while a light GAl 
is more consistent with HiRes-MIA in all the energy interval and with Auger 
between 4 — 6 x lO^'' eV. Between ~ 6 x 10^''' and ~ 2 x 10^^ eV Auger data 
is more or less in-between both solutions. These results are qualitatively inde- 
pendent of the HIM. Even if the assumption of a pure chemical composition 
for GAl and GA2 is artificial, it is clear that the present discrepancies between 
the various experimental results are of astrophysical significance, since they 
have quite different implications with respect to the nature of most powerful 
Galactic accelerators. However, and surprisingly enough, the uncertainty as- 
sociated with the several HIM currently in use in the literature does not pose 
a too large qualitative problem from the astrophysical point of view; that is, 
unless our present understanding of the hadronic interaction processes is not, 
somehow, fundamentally wrong. 

In sections §3.3.11 and §3.3.21 we determined the composition, as a function of 
energy, of the extragalactic spectrum (EG) and the two additional Galactic 
components (GAl, GA2) that fits X^ax data of Auger and HiRes, under the 
simplifying assumption of a binary mixture of p and Fe for the cosmic ray 
flux. 

The main result is that, regardless of the experimental data set considered, 
the composition has to be mixed to some extent along all the spectrum. The 
assumed HIM plays a very important quantitative role in the inferred compo- 
sition of any individual component. 

An important point to note regarding GAl and GA2 is that, even if they 
have mixed composition, the composition profile inside each one of them is 
similar in the sense that the individual fluxes are lighter at lower energies and 
then become heavier as the energy increases. This is a systematic effect that 
is quantitatively, but not qualitatively affected by the HIM used. Profiles such 
of those we obtained for GAl and GA2 are similar to what can be expected 
from different populations of SNRs immersed in differing environments. If 
this is correct, then probably only SNRs are required in order to explain the 
main part of the Galactic flux up to energies ^ 3 x 10^^ eV. In any case, 
GA2 is lighter than GAl, which could be an evidence for the existence of a 
minor contribution coming from a different source like, for example, inductors 
associated with compact objects. 

The estimated average atomic weight (see. Fig. [20]) varies widely depending on 
the assumed HIM. However, it can be seen that, for a given HIM, HiRes always 
implies heavier compositions on the Galactic side of the transition region than 
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Auger. Despite this systematic difference, both data sets show the same trend 
of a decreasing value of (A) as a function of energy. 

The opposite happens at higher energies on the extragalactic side of the flux, 
where Auger implies an increasing (A) as a function of energy regardless of 
the HIM, and HiRes predicts either a rather constant or a much more slowly 
changing composition. Also, for the latter data set, the uncertainties intro- 
duced by the HIM are qualitatively more important, since the high energy 
slope of the profile can change signs depending on the adopted interaction 
model. 

It must be noted that, for both data sets, there seems to be a break in the 
average energy composition profile at an energy around ~ 10^^ eV or slightly 
higher, which may further highlight the physical association between the ankle 
feature in the cosmic ray spectrum an the Galactic-extragalactic flux transi- 
tion. 



4 Conclusions 



We have analyzed the matching conditions of the Galactic and extragalactic 
components of cosmic rays along the second knee and the ankle. 

We have calculated the diffusive Galactic spectrum from regular SNRs using 
the numerical diffusive propagation code GALPROP. We used this first as- 
sessment of the Galactic cosmic ray flux to analyze the matching conditions 
with two alternative models for the extragalactic flux: a pure proton model 
[2j and a mixed composition model [3]. 

The first step was the matching of the observed energy spectrum by combining 
the Galactic and extragalactic fluxes. From this process, it becomes clear that 
additional Galactic components must be at play in reality. The minimum 
amount of additional components that allows us to satisfactorily reproduce the 
spectrum is either one or two. The pure proton model is the less expensive one 
with only one required additional component, GAl. The mixed composition 
model on the other hand requires, besides GAl, another components at still 
higher energies, GA2. 

If no other information is used, both theoretical models are indistinguish- 
able from the experimental point of view. Therefore, we also analyzed the 
effect of incorporating composition information, in the form of elongation rate 
data {{Xmax))- For this study, the shape of each one of the components is 
determined from the matching of the observed total energy spectrum. The 
elongation rate of the combined fluxes is then fitted to the HiRes and Auger 
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data sets by changing appropriately the composition as a function of energy 
of each one of the components, Galactic and extragalactic. In order to make 
the analysis simpler, only a binary mixture of p and Fe is considered. The 
main result is that the additional Galactic components, GAl and GA2, must 
have a mixed composition. Furthermore, inside each one of this components 
there is a progressive evolution of the composition from lighter to heavier as 
the energy increases. This is consistent with this components being originated 
in different populations of SNRs. Additionally, GA2 is globally lighter than 
GAl, which could indicate the possible existence of a minor contribution from 
another acceleration mechanisms, without a rigidity cut-off, operating at the 
highest energies. 

The uncertainties introduced by the HIM are always important from the quan- 
titative point of view. From the qualitative point of view, however, there are 
some results that seem rather independent of the HIM, like the diminution in 
average atomic mass along the low energy branch of the ankle and the exis- 
tence of a discontinuity in the slope of the energy profile of (A) around the 
mid ankle. These results seem also to be supported by both. Auger and HiRes. 

Our present results may be certainly considered as preliminary due to ex- 
perimental uncertainties and simplifications in the numerical approach. First, 
there is the paucity of data involved in the determination of the energy spec- 
trum in the region encompassing the second knee and the ankle and the diver- 
gence between (Xmax) measurements by different experiments below 3 x 10^^ 
eV. A proper experimental characterization of this very important region will 
likely have to wait until the release of the KAS CADE- Grande [77| and Auger 
enhancement data [72]. Second, there are arguable simplifying assumptions 
related with our diffusive treatment of the Galactic component at the highest 
energies which is, very likely, undergoing a change in propagation to a full 
ballistic regime. This can be somehow mitigated by the fact that Fe nuclei 
should still be diffusive inside this energy interval, while protons would only 
deviate importantly from the diffusive approximation at the highest energies 
considered for the Galactic flux. In any case, the importance of the transition 
region as a play ground for disentangling the Galactic and extragalactic cos- 
mic ray fluxes is unquestionable and considerable effort should be invested in 
its full experimental characterization and theoretical modeling. 
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